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ABSTRACT 

Sonic discharge coefficients are presented for two different 
geometry flow nozzles using nitrogen gas at high pressures (100 atm 

c p 

(lOOxlO'H/m )) where real-gas correctons are significant. Throat 
Reynolds number range extended up to 8x10^. Experimentally obtained 
coefficients for a nozzle with a continuous and finite radius of 
curvature agreed with those obtained analytically to within 0.2 
percent. Experimental coefficients for » long -radius ASME nozzle 
agreed to within l/4 percent to an empirical equation representing 
the most probable subscnic discharge CJefficient, 

INTRODUCTION 

The sonic flow nozzle is now being recognized as an accurate , 
precise , and eaoy to use head -type of flowmeter for controlling 
and measuring the mass flow rate of gases over a wide range of flow. 
References [1 - 5] 1 have pointed out the attractive and distinguishing 
features that make the sonic flow nozzle practical. Mass flow rate is 
easy to measure and control since it is directly proportional to up- 
stream pressure level. Upstream pressure tap error is negligible 
since the approach velocity is usually small. Absence of a differential 
pressure measurement eliminates the throat tap and simplifies construction. 
Accuracy and ease of use are superior to a subsonic nozzle because 
the hard to make measurement of the small differential pressure at high 
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pressures is not needed. Downstream pressure pulsations do not affect 
the mass flow rate because of sonic flow at the throat. 

The above advantages of a sonic nozzle are new being utilized 
[6 - 111 , However , none of the above references have reported experimental 
C D data at high-pressure levels where real -gas effects are significant 
(3 percent) and where the Reynolds numbers are as great as 1x10^, 

Such data would demonstrate the accuracy of the technique of combining 
the analytical values with the critical flow factors given by Johnson 
[12 - 15], 

The experiment described belcw involve two different geometry flow 
nozzles operated over a Reynolds number range of 5x10^ to QxlO^. These 
high pressure (100 atm (lOOxlO 1 5 N/m ) ) experiments add new data to the scant amowit 
reported in the literature cm sonic flow nozzles and expands the range 
of sonic Cp values into the regime where real -gas corrections are significant. 

An additional objective was to analytically calculate the C D for one of the 
above nozzles and compare the results to the experimental data, 

APPARATUS AM) TEST PROCEDURE 
Sctiic-Noz 2 le Geometries Selected 

The first flow-nozzle geometry selected (fig, l) had the traditionally 
recommended quadrant of an ellipse for the convergent section [l6]. This 
section was tangent to a straight throat section of length 0,6 times the 
threat diameter. Nozzle inlet -to -throat area ratio was 5.5. The second 

flow-nozzle geometry represents a class of wall profiles having a finite 
and continuous radius of curvature from the inlet through the throat 
plane. In particular, circular-arc profiles which have a constant 
radius of curvature and are a subset of this class have been examined 
in references [3, 7, and 17], Such profiles produce a ccxitinuous 
acceleration of the fluid along the wall which makes the analyti- 
cal calculat on of the sonic C D tractable and less inaccurate. 

The design radius of curvature of the nozzle shown by figure 2 was 
1+ times the throat radius (normalized radius of curvature, Rjj = *0. 

However, due to machining inaccuracies the radius of curvature of the nozzle 
utilised (fig. 2) was found to continuously increase in value from the 
throat plane to the inlet. The throat plane for this nozzle had a 

value of 2.3 . Nozzle inlet -to -throat area was 7. Nozzle threat 
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diameter for both of the above nozzles was a nominal 0.25 inch#* 
(6.35xlO"^m). Throat diameter to upstream pipe diameter ratio. 0. 
equaled 0.15 for both nozzles. And both nozzles were fabricated from 
free -machining itainleaa steel vltn a epecified maximum surface 

r J7 

roughness of 3«2xlO inches (8x10 m). 

Nitrogen-Gas Flow Facility 

A schematic drawing of the test facility is shown in figure 3. 
Nitrogen gas supplied from a portable trailer was throttled from 160 

C ry 

atm (l60xiCr newtons per square meter (N/;. 1- )) to the desired operating 
pressure levels which were varied fran 5 atm (5xl0^N/m 2 ) to 100 atm 
(lOOxlO^N/m 2 ) . Maximum flow rate was 3*71bm/sec (0.75 kg/sec). A 
bundle of tubes with each tube having a length-diameter ratio of 10 
was used just downstream of the pressure regulator to minimize any 
swirl in the flow, A sharp-edged orifice (modification of ref, 18) was 
used as the working standard with a maximum differential pressure of 
7 psi (0.48xl0 ? N/m ), The orifice, 60 -diameter -length approach pipe, 
and the 25 diameter -length downstream pipe were calibrated as an integral 
assembly on both water and air. Inaccuracy of the water stand used for 
this calibration was less than +0.2 percent. 

Differential pressure across the orifice was determined with a 
fused-quartz Bcurdon-tube precision pressure gauge [19], Biaccuracy 
of this gauge is less than +.05 percent of full scale. The null out- 
put from the electro -optical transducer of the gauge was displayed on 
an oscilloscope and time averaged. Differential pressure lines were 
matched for time delay find appropriate volumes and sintered porous - 
metal snubbers were introduced in each line to equalize time constants. 
This filter network reduced differential -pressure measurement errors 
due to absolute pressure -level fluctations which had a root -mean -square 
magnitude of approximately 9 percent of the differential pressure before 
filtering. Tbe finai orifice differential pressure inaccuracy was 
determined to be if . -c. then +0,1 percent. At absolute -pressure levels 

r— 0 

greater than 30 / 3 lx 10 N/tn - ) , the electro-optical transducer output 

of the gau*e was corrected by a method described in reference [20], 

The pipe section between the orifice and the nozzle test section 
as well as the nozzle test section itself were thermally insulated to 
give a uniform gas -temperature distribution across the nozzle inlet 
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plane. Orifice gas temperature was measured with a thermocouple prcfoe 
5.5 pipe diameters downstream of the orifice [21], Nozzle inlet gas 
temperature was measured with a thermocouple probe 3.*+ pipe diameters 
upstream of the nozzle inlet. Also, near this location the nozzle 
upstream static pressure was measured. Nozzle downstream pressure was 
measured sLout 1 pipe diameter downstream of the nozzle exit plane. 

The thermocouple electranotive force from the orifice and nozzle 
thermocouples were measured to a corresponding temperature inaccuracy 
of +O.L°F (+0.2K) . Operating absolute -pressure levels inaccuracies of 
less than +0.1 percent were achieved by using four different sets of 
gauges, of the Bourdon type, to span the pressure range used and by 
having the gauges periodically recalibrated. 

After instrument calibration and a slow pressurization of the facility 
to the maximum pressure level for the test to be run, system leak checks 
and instrument calibrations were performed. Nitrogen gas flow was then 
initiated by fully opening the throttle valve. The measured upstream to 
downstream pressure ratio across the test nozzle was always greater than 
2 which insured that the flow through the nozzle was indeed sonic. Not 
until the temperature of the facility had stabilized, approximately 5 minutes 
later, were data taken. Each data point consisted of a 1 minute run at 
constant flow rate. The flow rate was then changed by reducing the nozzle 
upstream pressure with the pressure regulator. 

ACCURACY AND REPEATABILITY 
Accuracy 

The overall relative inaccuracy associated with the value of C D is 
found by combining the relative uncertainties previously stated for each 
variable. Variables considered are the orifice flow coefficient, orifice 
and nozzle areas (and the effect of thermal contraction on throat area) , 
pressure levels, temperature levels, and the orifice differential pressure. 
Combining the relative uncertainties of the above mentioned variables 
for a throat Reynolds number of 1x10^ (upstream pressure equals 11 atm 
(lbcJ.O^N/m 2 )) yields on overall root -sum -square error of +.25 percent [22], 

Real -gas factor - This factor [12] corrects the conventional ane- 
dimensional, isentiopic flow relation used for calculating the main- 
stream mass flow rate by accounting for the variation of the ratio of 
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specific heat and the compressibility factor. Neglecting to use this 

critical flow-factor can give an error in mass flow rate exceeding 3 

5 2 

percent at a pressure level of 100 atm (100x10 N/m) for nitrogen gas 
[23]. 

An approximate value for the uncertainty in the critical -flow factor 
is obtained by comparing two different sources for nitrogen gas [12, 15]. 

The comparison shows at a temperature of 35°F (275K) and a pressure of 
10 atm (lOxlO'N/m ) the critical-flow factor from reference 12 exceeds 
the value from reference 15 by .03 percent. Whereas at 100 atm (100xKrN/m ) , 
the difference is 0,1+ percent. This difference increases monotonies Uy 
with pressure level. Accounting for this critical -flow factor uncertainty 
at a throat Reynolds number of 8x10^ (upstream pressure equals 95 atm 
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(95xKrN/m )) increases the overall root -sum -square relative error to 
+0.1+1+ percent. The critical -flow factors from reference 12 are used in 
this report. 

Repeatability 

A calibration consists of a set of data points ranging from the 
maximum to the minimum pressure level of any one of the four pressure - 
level ranges examined. Figure 1+ represents the total number of data 
points from 20 individual calibrations for the long -radius flow nozzle 
and figure 5 with 21 individual calibrations for the continuous wall 
curvature nozzle. At a Reynolds number greater than 1,6x10 . the 
repeatability is +0.2 percent for a 95 -percent confidence band for a 
single calibration [24], The confidence band for a single calibration 
is the band within which it is predicted that all the data points of the 
next calibration will lie with a 95 -percent probability. Stated band 
widths for a single calibration and for the mean are determined from 
t -distribution statistical calculations using the finite number of data 
points existing at selected Reynolds numbers. 

SONIC NOZZLE PERFORMANCE 
Experimental Results 

Accurate interpretation of sonic nozzle performance as revealed by 
the Cp curve shape and its level requires knowledge of the upstream 
velocity profile. As shown by Ferron [25]. the more uniform the upstream 
velocity profile the lower the C^ value. Inlet flow to both of the 
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nozzi.ee tested is regarded as fully "developed turbulent flow due to the 
nozzle upstream-pipe length -to -diameter ratio of 60. Also, the 
measured inlet -temperature profile variation is less than 2°F (IK) at 
the lowest mass flow rate. This essentially uniform temperature profile 
is due to the thermal insulation placed on the upstream pipe wall. 

Initial air calibration - Sonic flow nozzle C D variation with 
throat Reynolds number at flow conditions where leal-gas effects are 
small was determined for air by using the Lewis flow standards facility. 

The results are shown in figure s 4 and 5. The Reynolds number is based 
on the viscosity of the gas upstream of the nozzle and nozzle throat 
diameter. As seen from figure 4 , at the lcrw Reynolds numbers , the 
level of the curve agrees with the independent calibration 

obtained from the Lewis flow standards facility. This indicates 
systematic errors frcm instrumentation are negligible. 

Long -radius AS ME f low nozzle - A sufficiently high number of data 
points have been obtained that permits reliable calculation of the sonic 
nozzle C D mean curve shape and the 95 -percent confidence band for a 
single calibration and for the mean curve. The mean curve for the sonic 

C variation with Reynolds number for this nozzle is shown by figure 4. 

D 5 

It ranges in value from .983 at a Reynolds number of 3x10 to ,992 at 

a Reynolds number of 8x10^. Several features are indicated by the 

variation in the mean CL curve. First, starting at the low Reynolds 
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number, the C D value increases up to a Reynolds number of 1x10 which 

is most probably due to a decreasing boundary -layer thickness. Immediately 

following is a region of apparent transition from a laminar to a turbulent 

boundary layer. This region extends over a Reynolds number range of 
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3xl0 ? to 2xL0 . As can be seen from the increased size (+0.3 percent) 
of the 95 -percent confidence band for a single calibration, the C D value 
is relatively ill defined over this Reynolds number range. 

The apparent transition is most likely promoted by the local 
adverse pressure gradient generated by the discontinuous wall-radius of 
curvature at the Juncture of the convergent section with the cylindrical 
throat section [26], After this region of apparent transition, the C D 
curve monotonically increases in value with Reynolds number. At these 
Reynolds numbers the 95-percent confidence band for a single calibration 
is less than +0,15 percent. Also, the 95-percent confidence band for 
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the mean curve ia +0.05 percent. 

A comparison of the sonic nozzle C D mean curve to the beat empirical 
curve fit equation [27] for subsonic C D values is also shown by figure 4, 
Fran this comparison, the extent of valid representation, if any, for 
the sonic C D mean curve by this empirical equation is obtained. The 
empirical equation yields the most probable subsonic values for the 
long-radius ASME flow nozzle. This equation is limited to a maximum 
throat Reynolds number of 6x10^ at which the equals 0,993, Even 
though there exists a lack of similarity between the sub son! c c Q curve 
shape and the sonic mean curve shape, the most probable 

subsonic Cp values do represent the values to an inaccuracy 

of 1/4 percent for the experiments reported herein. 

Continuous curvature nozzle - the mean curve for the sonic 
variation with throat Reynolds number is shown by figure 5. It ranges 
in value from 0.989 at a Reynolds number of 6x10^ to 0,991 at a Reynolds 
number of 3x10^, A feature shown by the shape of the mean curve is 
the apparent transition for the Reynolds number region below lxKr*. 
Nevertheless , the exact shape of the mean curve is uncertain at these 
low Reynolds numbers where there are few data points. This 

uncertainty is indicated by the value of the 95 -percent confidence band 
for the mean curve increasing from +0,04 percent at a Reynolds 

* C 

number of 1x10 to +0.4 percent at a Reynolds number of 6x10 . For 

Reynolds numbers greater than 1x10^, the mean C D curve monotonically 
increases in value at a greater rate than the long -radius flew nozzle 
mean curve. Consequently, from the above mentioned factors for thiB 
sonic nozzle geometry, the flow range examined which is free of ill- 
defined transition regions is noticeably greater than the flow range 
for the long -radius flow nozzle. At these higher Reynolds numbers, 
the 95 -percent confidence band for a single calibration is less than 
+0.2 percent. 


Comparison to Analytical Results 

Nozzle flow model - Only sonic nozzles with continuous and finite 
radius of curvatures are considered herein. The gas flow through the 
nozzle is divided into two different flow regions. One is the mainstream 
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gas -flow region that Is treated as an irrotat iomal , isentropic flow of a 
perfect gas. Hi this region, fluid acceleration near the wall and the 
resulting inertial forces produce a nan -uniform velocity profile 
(sonic line curvature) at the throat. Vincent [173 shows that inertial 
effects can affect the C^ by as ouch as 1 percent. The second is the 
viscous flow region along the wall occupied by the boundary layer. 

The boundary layer can lower the value by as much as 5 percent 
at low Reynolds numbers. The effect of aach of these two regions on 
the value of the sonic is separately analysed. 

The mainstream gas flow parameters across the nozzle throat plane 
are calculated by using q computer program based upon reference 28 for 
axisymmetric flow. Calculation results are dependent upon the following 
input variables: (l) normalized wall radius of curvature at the throat 

(Rjj) ; ( 2 ) nozzle throat radius; (3) nozzle axial length; (4) specific - 
heat ratio of the gas; ( 5 ) boundary -layer displacement thickness at the 
throat. The normalized throat radius of curvature of the nozzle studied 
herein is equal to 2,3 , Actual wall profile, which is not a circular 
arc, is piecewise curve fitted by the program from reference C? 93 which 
enables the wall Mach number distribution to be determined. 

Two-dimensional , boundary -layer parameter values are determined 
by a computer program from reference [3C], The above mentioned wall 
Mach number distribution serves as part of the input data to this 
boundary -layer program. Laminar calculations are based upon the 
work of Cohen and Roshotko [ 3X3 . And the work of Sasman-Cresci [ 32 ) 
is used as the basis for the turbulent boundary -layer calculations. 
Modifications to this program with real -gas subprograms of reference [ 33 ] 
permit the accurate calculation of due speed of sound, gas density, 
and specific heat of nitrogen gas. Furthermore, transport property variation 
with pressure level is also accounted for by use of data from references 
[34-36]. The resulting boundary -layer displacemnt thickness at the 
throat is then used as part of the input to the above isentropic gas- 
core flew program for calculating the final value of the acnic C^. 

Analytical results comparison - The analytical sonic C^ variation 
with throat Reynolds number for the continuous curvature nozzle is shown 
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by figure 6 for nitrogen gas for both laminar and turbulent boundary- 
layer conditions. Also shown is the mean sonic C D curve for the real- 
gau nitrogen flow datM, Agreement between the analytical and experimental 
curves over the entire Reynolds number range experimentally investigated 
is good considering the previously discussed uncertainty in the mean 
Cp curve at the low Reynolds numbers. However, the analytical 
variation for the laminar boundary layer differs from the Lewis in- 
house flow facility calibration by as much as 0.6 percent at the low 
Reynolds numbers. The difference between the analytical curve for the 
turbulent boundary layer and the mean curve is equal to or less than 
0,2 percent. This difference is within the mean C D curve 95 -percent 
confidence band for a single calibration, A comparison of the level 
and the slope of the analytical curves to the mean C Q curve tends to 
verify that the boundary layer is turbulent over the greater extent of 
the Reynolds number range covered. As previously stated, the 
experimental value is determined by using a^critical flow factor. 
Considering the previously discussed sources of error, it is seen that 
the sonic for the nozzle tested csn be determined to practical 

flow metering accuracy by the results from the turbulent boundary -layer 
calculations including real -gas effects. Therefore, the sonic for any 
nozzle with a continuous and finite radius of curvature (R^ >2) from 
the inlet to the throat plane can be expected to be represented over 
its practical Reynolds number range by the above analytical method 
for calculating the turbulent boundary layer. This assumes a corner 
exists at the inlet to augment early transition and thereby extend 
the useful Reynolds number range. Also, the installation of the 
sonic nozzle should be according to recognized standards. The 
ability to accurately calculate the is highly desirable because 
it provides the sonic flow nozzle with the potential for becoming 
a compressible flow standard. 
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Figure 3. - Nitrogen gas flow facility tor sonic (low nozzles. 



Figure 4 - Experimental discharge coefficients tor the tong-radius 
ASME (low nozzle using nitrogen gas. 
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Figure 5. - Experimental discharge coefficients for the continuous 
wall-curvature sonic flow nozzle using nitrogen gas. 



Figure 6. - Comparison of the analytical results to the experimental 
mean-value curve for the continuous w3ll-curvature sonic flow 
nozzle. 
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